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Chromatin looping is a major epigenetic regulatory mechanism in higher eukaryotes. Besides its role in 
transcriptional regulation, chromatin loops have been proposed to play a pivotal role in the segregation of 
entire chromosomes. The detailed topological and entropic forces between loops still remain elusive. Here, we 
quantitatively determine the potential of mean force between the centers of mass of two ring polymers, i.e. 
loops. We find that the transition from a linear to a ring polymer induces a strong increase in the entropic 
repulsion between these two polymers. On top, topological interactions such as the non-catenation constraint 
further reduce the number of accessible conformations of close- by ring polymers by about 50%, resulting in 
an additional effective repulsion. Furthermore, the transition from linear to ring polymers displays changes 
in the conformational and structural properties of the system. In fact, ring polymers adopt a markedly more 
ordered and aligned state than linear ones. The forces and accompanying changes in shape and alignment 
between ring polymers suggest an important regulatory function of such a topology in biopolymers. We 
conjecture that dynamic loop formation in chromatin might act as a versatile control mechanism regulating 
and maintaining different local states of compaction and order. 



I. INTRODUCTION 

Ring polymers are abundantly found in nature, com- 
mon examples range from the DNA of bacteria to certain 
virus types''^, mitochondrial DNA in higher eukaryotes 
and proteins^ . These findings suggest a structural advan- 
tage of ring polymers over linear ones. The genomic con- 
tent of higher eukaryotes inside the interphase nucleus, 
however, has a much more complex organization. Com- 
paction is accomplished on several scales, the detailed 
fol ding pathways being closely related to genome func- 
tiorPSl. Nowadays, abundant experimental and theoreti- 
cal evidence suggests that chromatin looping acts as an 
important link between chromatin organization and func- 
tion. The existence of chromatin loops on scales from a 
few kilo basepairs (kb) up to several mega basepairs (Mb) 
has been evaluated by3C/4C/5C experiments!-^-'-'. Several 
mechanisms for looping have been proposed: Amongst 
others, it was suggested that these loops arise when dif- 
ferent genes and regulatory elements co-localize in tran- 
scription hubs or factories^, where a high concentration 
of polymerases engages transcription. 

Inconsistent with theoretical predictions for solutions 
of linear polymers, chromosomes arrange such that they 
occupy distinct territories, the overlap volume being 
quite small. A similar effect has been found even for 
small parts of single chromosomes spanning several mega 
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basepairs (MbjP. In principle, active ATP-consuming 
mechanisms might be responsible for maintaining segre- 
gation. Rosa and Everaers^*' proposed that chromosomes 
are aggregating into distinct compartments because of 
their large relaxation time, making it impossible for them 
to intermingle during the time of one cell cycle. How- 
ever, this study does not explain segregation of smaller 
regions within chromosomes. Vettorel and co-workerJSl 
suggested that the collapsed state might be responsible 
for segregation: Polymer gels, where single polymers are 
in a collapsed state, show no reptation, therefore inter- 
mingling between different polymers becomes practically 
impossible. Recently, evidence from polymer models has 
indicated that chromatin loops force such a compar tmen- 
talization already by virtue of purely entropic force^^^*^. 

While the mean squared radius of gyration of iso- 
lated ring polymers and a system of two catenated ring sSl 
displays a scaling similar to that of a self-avoiding walk, 
where (i?^) ^ N'^^ with v « 0.58#^, non-catenated rings 
in a melt become compact adopting a scaling exponent 
of = l/ ji^l^'^ l. In fact, chromatin organization is not 
consistent with a simple melt of ring polymers. Scaling 
exponents found in experiments differ a lot, ranging from 
1/ « 0.1 - O.^Slto V = Q.^, the level of compaction de- 
pending on gene activitjESl_ Dynamic chromatin loops, 
however, have been shown recently to provide a consis- 
tent framework of chromatin organizatiorP^. 

The purpose of this study is to deepen the understand- 
ing of the topological interactions ring polymers exert 
on each other as well as the changes in their conforma- 
tional properties induced by such forces. While topo- 
logical effects of rings have been considered in several 
studie^ ^ l ^^ l ^^ l, the induced forces have not been quanti- 
tatively assessed. A system of two ring polymers can be 
viewed as a toy system for the influence of loop forma- 
tion on chromatin folding. While the biological system is 
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beyond doubt much more comphcated - multiple loops 
and the dense system of polymers playing a dominant 
role - the effect of topological forces and the advantage 
of the ring shape for biopolymers can be highlighted best 
by studying it apart from other influences. We evaluate 
the potential of mean force between two ring polymers 
dependent on their mutual topological state to derive a 
quantitative measure for the entropic repulsion two rings 
or loops exert on each other. Therefore, Monte Carlo sim- 
ulations of ring polymers with lengths of up to TV = 2048 
monomers are conducted. We separate the entropic ef- 
fects due to excluded volume interactions from the topo- 
logical interactions arising from the non-catenation con- 
straint. The ring polymers' conformational properties 
and their mutual alignment subjected to the topological 
and entropic forces are investigated. In fact, we find that 
ring polymers display a markedly stronger repulsion and 
adopt a more ordered state compared to linear polymers. 
The findings indicate a natural benefit of rings over lin- 
ear ones, leading to the conjecture that chromatin loops 
not only facilitate entropy-driven segregation of chromo- 
somes or intra-chromosomal regions, but act as a versatile 
control mechanism regulating and maintaining different 
states of compaction and order. 

The paper is organized as follows: In section |ll] we 
revisit the theory behind the effective potential, present 
details of the simulational method and the numerical cal- 
culations of the effective potential. In section III A 



we 



calculate the potential of mean force for linear polymers 
and review scaling theories. We turn to ring polymers 
in section |IIIB| and study the effect of topological inter- 

The following two sections are 



actions in section IIIC 



dedicated to the investigation of how the dimensions and 
shape of a ring changes when bringing two rings into 
contact. Finally, in section |IIIF[ we study the mutual 
alignment of rings considering the effect of topological 
constraints. 



II. SIMULATIONS AND METHODS 
A. Effective potentials 

To investigate the potential of mean force exerted be- 
tween the centers of mass of ring polymers, we use the 
concept of the effective potential. The method has been 
used for a variety of polymer systems, and f or som e even 
analytical approximations have been found I^SIlli. A rec- 
ommendable review is given in Ref.l^. 

Let the monomer positions and momenta of the two 
polymers be denoted by rf and pf (i — 1, . . . ,N,a — 
1,2), respectively, a denotes the polymer chain index 
and i indexes subsequently the monomers of the polymer. 
The properties of the system are given by the complete 
Hamiltonian 
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For the problem of ring polymers investigated here, the 
interaction potential V^({rf }) is given by an excluded vol- 
ume term and eventually by a term describing the topo- 
logical constraints. 

The potential of mean force can be gained via the par- 
tition sum of the system constrained to a fixed center-of- 
mass, 
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The trace denotes the integral over all momenta and 
monomer coordinates. 

The effective potential is now defined by the relation 

exp(-/3i7eff(R)) = Ze(R) 

For purpo ses of calculation, however, the following rep- 
resentatiorP^I^ of the effective potential is more useful 



t/eff(i?) =-fcBTln 



ZcjR) 



(1) 



where Z\ is the partition sum of a single polymer. Equiv- 
alently, Z\ may be imagined as the partiton sum of a 
two-chain system whose centers of mass are infinitely far 
apart. 



B. Monte-Carlo algorithm 

The generation of ring polymers is accomplished using 
a well-tested latt ice M onte Carlo algorithm, the bond 
fluctuation modepS122l. Simulations are set up on a cu- 
bic lattice, where each monomer of the ring resides in 
the center of a unit cube and occupies all its eight ver- 
tices. Each monomer is connected by a bond to its two 
neighbors along the contour of the polymer, ensuring the 
connectivity of the chain. The lattice constant is taken 
to be unity and all measures of length in the following 
are given in multiples of the lattice constant. The length 
of the bonds between two neighboring monomers are al- 
lowed to fluctuate, but constrained to a set B of 108 
vector^2ll One Monte Carlo trial move encompasses ran- 
domly selecting one monomer and proposing a move to 
one of its nearest neighbors on the cubic lattice. Using 
these local moves combined with the constrained set of 
bond vectors it has been shown that no bond-crossings 
can occur during the simulation ruiP^, thus the topolog- 
ical state of the rings is preserved while propagating the 
system. Excluded volume interactions are taken into ac- 
count by preventing a lattice site to be occupied by more 
than one monomer. 

Isolated ring polymers with chain lengths between 
= 32 and N — 2048 are simulated. As only local moves 
are applied, consecutive conformations are highly corre- 
lated. To determine the relaxation time of a ring poly- 
mer, we evaluated the autocorrelation function C{t)'^^ of 
the squared radius of gyration A{t) = R'i(t) and from 
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a center-of-mass distance r is then given by 



FIG. 1. Illustration of different topological states of two ring 
polymers. The left figure shows a non-catenated conformation 
with Gauss linking number 4> = 0. On the right, a catenated 
conformation is displayed having a linking number of "I> = 1. 



C{t) we estimated the integrated autoc orrelation time r 
using Sokal's windowing procedur^SEH^ Two subsequent 
conformations are considered independent after 5r Monte 
Carlo steps (MCS). In the following we define one MCS 
as N trial moves, where N is the chain length and there- 
fore the total number of monomers in the system. Thus, 
in one MCS each monomer is proposed on average one 
move. For each chain length at least 10 000 independent 
conformations were produced. 

Simulations are performed on a lattice with periodic 
boundary conditions, but the algorithm keeps track of 
unfolded coordinates, such that the entire configuration 
space can be explored. The linear lattice size L is cho- 

1/2 

sen large enough, i.e. L ^ (^g) ; so that monomers 
do not come into contact in the lattice space by peri- 
odic backfolding of coordinates, otherwise biased confor- 
mations would result. Lattice sizes used vary between 
L = 256 up to L = 800, depending on chain length N. 



C. Calculation of the effective potential 

Single ring conformations are used to calculate the av- 
erage interaction energy between two polymers by means 
of equation (Tj). This procedure is we ll-established and 
has been applied in several studie^^HUi] j^^ first, the set of 
sampled conformations C is split into two parts of equal 
size. Then, we randomly pick two conformations, one 
from each subset. These two ring conformations are then 
positioned such that their centers of mass have a certain 
distance r, the angular positioning in doing so is cho- 
sen randomly. In order to ensure all monomers residing 
on lattice sites after shifting, a maximum distance in- 
accuracy of Sr « 0.87 lattice units is accepted. In the 
next step, the interaction energy Ei{r) of the two-ring 
conformation and the corresponding statistical weight 
Wi{r) = exp{—Ei/kBT) are calculated. 

Repeating this operation for a number K of randomly 
selected pairs of conformations, the effective potential at 
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The standard error of the effective potential is deter- 
mined by randomly subdividing the K two-chain con- 
formations into AI = 50 smaller subsets and calculating 
the effective potential C/m(r) for these subsets^^-. The 
standard error is then calculated by 



AC/„a(r) 
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In a first step, only taking into account excluded vol- 
ume interactions, the statistical weight Wi{r) of a two- 
chain conformation is set to Wi = in case the excluded 
volume condition is violated, i.e. two monomers occupy 
the same lattice site, and Wi = 1 in case it is not. 

However, care has to be taken using the bond fiuc- 
tuation model. During the simulation run it is ensured 
that no bond crossings can occur; this, however, is not 
guaranteed if we just shift two conformations inside each 
other. There are two possible combinations of bond vec- 
tors, where the excluded volume condition is not vio- 
lated, but where the bonds come into contact. Consider 
the two bonds spanned by the following four monomers 
{(0,0,0) ^ (3, 1,0), (1,2,0) ^ (2,-1,0)}. Both bond 
vectors (3, 1, 0) G S and (1, —3, 0) £ B are valid and the 
excluded volume condition is satisfied. The second prob- 
lematic case is {(0, 0, 1) ^ (2, 2, 0), (0, 2, 0) ^ (2, 0, 1)}. 
As can be easily seen, the participating beads are stuck 
in these positions, there is no valid move for any of the 
beads using the restrictions of bond vectors the algorithm 
is subjected to. Therefore, during a simulation run such 
a situation can never happen, assuming the start config- 
uration is chosen properly. Here, we test each two-chain 
conformation, after positioning them with their centers 
of mass a given distance r apart, on such a situation. If 
a bond crossing occurs, a weight of Wi = is assigned to 
this conformation. 

In this study, we are furthermore interested in the 
topological state of a two-ring conformation. We cal- 
culate the weight factor differently whether we consider 
a preserved topological state or not. In the case that we 
force a certain topological state, we have to add another 
criterion for accepting a two-chain conformation (i.e. as- 
signing a weight Wi ^0). Two rings that are positioned 
at a centcr-of-mass distance r might be non-catenated or 
catenated. Furthermore, the degree of catenation might 
vary (see illustrations in Fig. [T]). We determine the topo- 
logical state by a topological invariant: the "Gauss link- 
ing number"—. This invariant has been used in analyti- 
cal stud ies on interlinked rings to keep track of a certain 
topologjESED. It is defined by 
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The closed line integrals are evaluated along the contours 
of the two rings, denoted by Ck {k = 1,2). The vector 
function = J^k{s) denotes the three-dimensional coor- 
dinates of the ring polymers, parameterized e.g. by the 
contour length s. 

The Gauss link invariant has the advantage that it is 
easy to evaluate, although it has the disadvantage that it 
is not one-to-one. The "Whitehead link"'"9 is one exam- 
ple for a two-ring conformation having the same linking 
number as non-catenated rings. However, for this link, 
self-intersections of one ring are necessary, which do not 
arise during the simulation runs. Therefore we consider 
the Gauss linking number as an excellent quantity to de- 
termine the topological state. The result of the integral 
is an integer for non-overlapping chains, especially $ = 
for non-catenated chains and $ = 1 for simple catenated 
chains (Fig. [ij . We numerically evaluate this integral to 
classify the topological state of a two-ring conformation 
and assign a weight Wi = for all two-ring conformation 
not in the topological state we are interested in. Thereby 
we are able to calculate the effective potential separately 
for any given topological state. 

III. RESULTS 

A. The effective potential of self-avoiding walks 

We first recapitulate the results for the effective po- 
tential between the centers of mass of two self-avoiding 
walks both to validate the algorithm used and to relate 
the strength and form of the potential for linear polymers 
to their ring counterparts. 

A simple mean field argument can be devised for the 
entropic repulsion of two polymer coils with excluded vol- 
ume at full overlap, i.e. with zero center-of-mass sepa- 
ration. It has been first proposed by Flory and Krig- 
baunPH and later been corrected in other studieJ^SEH. 
The overlap volume is given hy V ^ Rg where Rg is 
the radius of gyration of one polymer coil. A mean-field 
argument yields the effective potential at full overlap^. 
Scaling argument^^, which take into account the corre- 
lation hole effeciP^, show that the probability pi of one 
monomer of chain A being in contact with any monomer 
of B is given by pi = (na^) ^"-^ . Thus the probability 
that none of the monomers are overlapping is given by 
p = (1— pi)^ = (1— j^)^ ■ The effective potential can 
then be calculated as the logarithm of this probability, 

U{r = Q)^-N\n{l-^-^ (3) 

In the limit of infinite chain length iV oo, the effective 
potential at full overlap approaches a constant value in 
the order of fcsT. Consequently, linear polymer coils with 
excluded volume have a rather soft potential, allowing for 
a high level of mutual penetration. 

Results for the self-avoiding walk are displayed in 
Fig. [2] The effective potential ?7sAw(^) in Fig.[2]\ reveals 



a decrease with growing chain leiigth, inconsistent with 
the mean-field argument by Flory^ but consistent with 
eq. ([3]). The effective potential at full overlap is shown 
in a scaling plot (cf. Ref.l^ in Fig. [2|3. Data points are 
on a straight line as expected, deviations are most likely 
due to lattice effects: When the two polymer coils are 
brought to a certain distance, the nearest lattice site is 
chosen, such that the distance can deviate from the target 
distance r — Q. These deviations are more pronounced 
the smaller N . 



B. Rings have a stronger repulsion than linear polymers 

To study the effective potential between ring polymers, 
we first neglect topological constraints between the two 
rings when bringing them in proximity. However, each 
single ring obeys the topological constraint of unknot- 
tedness by virtue of the simulational method. Thus the 
effective potential is given by the logarithm of the proba- 
bility that the ring conformations do not occupy at least 
one lattice site in common. 

FigurejsjA shows the effective potential t/ring(?') for ring 
polymers where mutual topological constraints are ne- 
glected, i.e. the linking number between the two ring 
polymers is arbitrary. The effective potential again is 
dependent on chain length iV, however, differences be- 
come smaller for large chain lengths. The mean field 
argument for the effective potential at full overlap estab- 
lished in eq. ([3| does not depend on the connectivity of 
the chain. Indeed, it is not only valid for linear poly- 
mers, but also for rings. Figure [3j3 shows the effective 
potential at full overlap Luring (0) such that a linear ex- 
trapolation to — > 00 becomes possible. Clearly, the 
effective potential adopts a constant value in this limit. 

Interestingly, the repulsive interactions between two 
ring polymers are much stronger than for their lin- 
ear counterparts of equal chain length N . The ratio 
C^ring(?')/t^SAw(^) IS displayed in Fig.|4]and has a value of 
about 3 for center of mass (CM) separations below 2i?g. 

C. Topological constraints induce additional repulsion 

In the last section we have shown that the repulsive 
forces acting between two ring polymers whose centers of 
mass are brought in close proximity are much stronger 
than for linear self-avoiding walks, the potential being 
about 3 times larger. This rise in strength results from 
the increased density of monomers induced by the ring 
structure. Here we focus on the additional effect topo- 
logical constraints induce on the potential of mean force 
between the polymers. Much attention has been paid to 
the influence of topological constraints on ring polymers. 
The interaction between non-linked ring polymers with- 
out excluded volume has been investigated for chains of 
length up to = SCpS. In a recent studjEl we have 
analyzed the behavior of two catenated rings. It turned 
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FIG. 2. Effective potential acting between the centers of mass of two self-avoiding walks. A. This panel displays the effective 
potential UsAw{r) for different chain lengths. The a;-axis is scaled by the root mean squared radius of gyration Rg of isolated 
chains. Errors are smaller than the point size and therefore not shown. The effective potential is dependent on chain length. 
B. Scaling plot of the effective potential at full overlap. The effective potential L''saw(0) reaches a finite value for large chain 
length as expected from scaling theory (cf. equation ([3|). 




FIG. 3. Effective potential ?7ring(r) of ring polymers disregarding their mutual topological state. A. This panel displays the 
effective potential in units of ksT for several chain lengths. Errorbars are not shown as they are smaller than the point size. 
For reasons of comparison, the a;-axis is scaled with the radius of gyration of the corresponding isolated ring polymer. B. The 
scaling law for the effective potential at full overlap follows the theory from eq. Q allowing for a linear extrapolation to the 
infinite chain limit when plotting e~'^'°'''^ vs 1/N. For the extrapolation, only chain lengths larger than A'^ = 256 are used. 



out that while the dimensions of catenated rings show 
the same scaling behavior than their hnear counterparts, 
the shape changes significantly. The topological states 
of ci rcula r DNA have been analyzed in various publica- 
tiond^^. The studies on non-catenated and unknot- 
ted ring polymers in the m elt revea led that their di- 
mensions change dramaticalljffillfillll They behave like 
compact polymers with a radius of gyration scaling with 



N^^^. The dominant entropic force driving this com- 
paction stems from the topological constraint of non- 
catenation. To our knowledge, a quantitative study of 
this topological potential has not been considered yet. 
In the framework of effective interactions, the topologi- 
cal non-catenation constraint can be included easily. As 



explained in the section II C we determine the topologi- 
cal state of a two-ring conformation by means of its Gauss 
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FIG. 4. Comparison of the effective potentials for ring poly- 
mers and self- avoiding walks. The ratio Uringir) /UsAw{r) is 
plotted against the center-of-mass distance r. The a;-axis is 
scaled by the root mean squared radius of gyration Rg of iso- 
lated chains for reasons of comparison. 



linking number $ (eq. ([2j)). A Gauss linking number of 
zero indicates the non-catenation of the two ring poly- 
mers. To obtain the potential of mean force under the 
non-catenation constraint, we calculate the fraction p of 
two-ring conformations with both zero linking number 
and fulfilled excluded volume constraint. The effective 
potential is then given by U^^^^(r) = — Inp (the index 
indicates the Gauss linking number $ = 0). The re- 
sulting potential U^^^^^{r) is shown in Fig. [sj^. We find 
pronounced deviations from a Gaussian shape for small 
CM separations r < IRg- 

From the effective potential C/°ing(^) '^^ non-catenated 
rings and C/ring(?') of rings without topological con- 
straints, the topological potential J/topl?"), i-e. the part of 
^ring('') stemming purely from the non-catenation con- 
straint, can be obtained. 



V(r-) = C/°ng(r) - [/H„g(r) 



The topological potential is shown in Fig. [5j3. It has a 
maximum at about r/Rg 0.55, a value which basi- 
cally does not change with chain length A''. For larger 
CM separations Uto-pir) drops to zero, as less and less 
two-ring conformations will be catenated when increas- 
ing r. There is also a decrease in the topological po- 
tential when going to small CM separations r < O.bRg. 
This decrease becomes smaller for larger chain length 
N. The reason for this behavior is that for small r 
(and small iV), two- ring conformations not satisfying the 
non-catenation constraint, often do not satisfy the ex- 
cluded volume constraint either, such that the confor- 
mations are not counted for the purely topological in- 
teraction C/top- However, more interesting is the maxi- 
mum repulsion exerted by the non-catenation constraint. 
We find that the topological potential is asymptotically 
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FIG. 5. A. The effective potential (7°iiig(r) for non- 
catenated rings. B. The topological potential Utopir) for 
non-catenated rings is obtained as the difference between 
U^\ng{r) and L'ring(''). It displays the potential arising purely 
from the topological non-catenation constraint. The inset 
shows the maximum of the topological potential depending 
on the chain length A'^. The maxima are obtained by a fit of 
t''ring('') to a third-order polynomial in the interval [0.3 : 0.7]. 



(see inset of Fig. [5|3). In numbers, this means that only 
about 48% (exp(— 0.65) ~ 0.52) of all two-ring confor- 
mations with fixed center-of-mass distance which would 
be allowed with respect to excluded volume are actually 
allowed with respect to topological constraints. In other 
words, only half of the two-ring conformations satisfying 
the excluded volume constraint are actually accessible 
when the non-catenation constraint is considered. Thus, 
the non-catenation constraint decreases the number of 
accessible conformations at small CM separations con- 
siderably. 

Interestingly, the topological interaction is on the order 
of magnitude of 1 k^T. This value was used for a crude 
estimate of the free energy of rings in a melt by Gates and 
DeutsctP^, where they assumed that the free energy cost 
for a contact between two rings is about IfcsT, hence 
in a three-dimensional melt the topological contribution 
to the free energy is F w ksTR^ /N . Although highly 
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FIG. 6. A. Mean squared radius of gyration for self-avoiding 
walks (SAW), ring polymers without mutual topological in- 
teractions and non-catenated rings vs. distance between the 
centers of mass. Radii of gyration are scaled by the corre- 
sponding values Rg of the isolated chains. B. Scaling of 
the radius of gyration with chain length. The leading term 
N"^" of a self-avoiding walk (y — 0.588) is divided out to see 
deviations from the linear chain behavior. 



simplistic, this scaling exponent v = Ijh was later re- 
produced by simulational resulttP^ more recent results, 
however, suggest that this behavior is only a cross-over 
to the behavior of compact lattice animals^. 



D. Dimensions of two-ring conformations 

When two polymers are brought closer and closer to- 
gether, lowering the distance between the centers of mass, 
not only the repulsive interaction increases dramatically. 
Also the conformational properties of the polymers are 
affectecP in an effort of minimizing the free energy. Here 
we investigate how the dimensions of a two-ring confor- 
mation change compared to a self-avoiding walk when 
approaching isolated polymers. 

There are two measures of the dimensions of the chain, 
which, however, exhibit the same scaling: the mean 
squared radius of gyration and the mean squared 

ring diameter ^d^) (see Ref-^SI). Figure|6j'V shows the root 
mean squared radius of gyration Rg{r) = ^ {Rg{rY) for 
self-avoiding walks (SAW), ring polymers without mutual 
topological interactions (arbitrary linking number <i>) and 
non-catenated rings (linking number $ = 0) in depen- 
dence of the scaled center-of-mass separation s = r / Rg. 



FIG. 7. A. Ratio between the mean squared radii of gyration 
of rings without mutual topological interactions (NT) and 
non-catenated rings (NC) vs. scaled distance r/Rg between 
the centers of mass. Results are shown for a chain of length 
A'' = 1024. B. The ratio is plotted against the inverse chain 
length allowing for an extrapolation to the infinite chain limit 
A' — > oo. Data is shown for two CM separations: s = and 
s = 1 to analyze whether deviations remain in the limit of 
large chain length. 




1 1.5 2 2.5 3 
CM separation r/Rg 



FIG. 8. Ratio between the mean squared radius of gyration 
and mean squared ring diameter vs center-of-mass separation 
r. Data is shown for a chain length of A^ = 1024. Results are 
shown for two-ring conformations without mutual topological 
interactions (open circles) and non-catenated ring conforma- 
tions (solid red circles). 
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The data are scaled with the radius of gyration for iso- 
lated polymers Rg = Rg(r = oo) (i.e. the radius of gy- 
ration of a polymer at infinite distance away from the 
second one) for reasons of comparison. For small CM 
separations (r < Rg) the radius of gyration strongly 
increases and its dimensions are more than 10% larger 
than in the isolated case. This is consistent with the re- 
sults for self-avoiding walks by Dautenhahn and Hah'^'^. 
The increase in dimension at separations r <C Rg origi- 
nates from the need of the chains to create space for the 
monomers of the second chain. At intermediate separa- 
tions of about r/Rg « 1.5 — 2, the radius of gyration at- 
tains a smaller value, thus the polymer effectively under- 
goes compaction. Both the compaction at intermediate 
separations as well as the swelling at high overlap remain 
in the limit of infinite chain length. An extrapolation of 
the radius of gyration at full overlap {s — Q) to N co 
shows that Rg{r = 0)/Rg — )■ 1.109(1) for rings without 
topological interactions and Rg{r = 0)/Rg — >■ 1.111(1) 
for non-catenated rings (see Supplementary Figure lH2l). 

One might ask, whether this change in dimensionality 
is accompanied by a change in the scaling law. We know 
from several studies that isolated ring polymers follow a 
similar scaling law as linear polymers {R'^) ^ N'^'^ where 
v w O.SJP^lSlIlI] However, non-catenated rings in a melt 
behave like compact polymers with a scaling exponent 
z/ = l/ jii | i'i' l Figure [6p shows the scaling of the radius of 
gyration with chain length in a log-log plot for two center- 
of-mass separations: r = (full overlap, region where the 
chain is swollen) and r — 1.7 Rg (region where the chain 
is compacted). The leading order term A^^" of the self- 
avoiding walk behavior (i/ = 0.58^^ is divided out. The 
data adopt a constant value in the limit of large chain 
length A^ independent of CM separation, indicating that 
the scaling behavior is independent of separation. As this 
also applies to rings with the non-catenation constraint, 
the compaction found in a polymer melt of non-catenated 
rings seems to require several chains, exerting forces from 
all sides on the polymer coil. 

Our results show that differences between the radius of 
gyration {RI)^^ of rings without topological constraints 
and of rings obeying non-catenation are rather 

subtle. The ratio of the radii of gyration {RI)^^ / {^1}^^ 
with respect to the CM separation is shown in Fig. ^p. 
for a chain length of A^ = 1024. For small separations 
r < 0.5i?g non-catenated rings are smaller, while for in- 
termediate separations (r ~ l^g) they are larger. The 
dependence of these results on the chain length A^ is given 
in Fig. [7f3 for two different CM separations. The dif- 
ferences are asymptotically stable both in the regime of 
small as well as intermediate separations. At full overlap 
(s = 0), non-catenated rings are larger in size by about 
0.8% in the asymptotic limit N od, for intermediate 
separations (s = 1), rings without topological constraints 
are larger by about 1.5%. Similar results are found for 
the mean square ring diameter (data not shown). 

As isolated ring polymers as well as self-avoiding walks 
only have one length scale, different measures of dimen- 



sion should yield a constant ratio. For a random walk, 
the ratio between the mean squared radius of gyration 
and the end-to-end distances is 1/6. For isolated ring 
polymers, the ratio (Rg) / (cP) extrapolates to 0.3053(2) 
in the asymptoti c limi t, while catenated rings display a 
ratio of 0.2995(3)fSIlIl. For the case studied here, the 
ratio between the radius of gyration and ring diameter 
/ (^cP') changes significantly for different CM sepa- 
rations (see Fig. 8]). From its isolated chain value the 
ratio decreases when bringing the centers of mass closer 
together. Extrapolation to infinite chain length N oo 
shows that the differences remain even in this limit (Sup- 
plementary Figure At full overlap (s = 0) the 
ratio extrapolates to (i?^) / {(P) = 0.2874(3) for non- 
catenated rings and / {(P) — 0.2891(3) for rings 
without topological interactions. 



E. Shape of two-ring conformations 

A typical measure of the shape of a polymer coil is its 
gyration tensoi-^^*^. It is defined by 



1 



N 



^mn — / ^ ' rn ' n 



(4) 



Here r*^*^ is the coordinate vector of the i^^ monomer 
and the subindex denotes its cartesian components. The 
matrix S is symmetric and positive semi-definite, thus it 
can be transformed to a diagonal matrix where the three 
eigenvalues Ai < A2 < A3 give the squared lengths of 
the principal axes of gyration of the associated gyration 
ellipsoid. In fact, the sum of the eigenvalues gives the 
squared radius of gyration i?^ = Ai -I- A2 -f A3. 

In a recent studjff^ on catenated rings we have found 
that the shape of catenated rings differs significantly from 
the shape of isolated ring polymers. This shows up in the 
ratios of the averaged eigenvalues (A3) : (A2) : (Ai), which 
are obtained from extrapolation to the asymptotic limit 
as 8.23(2) : 3.22(2) : 1 for isolated rings and 8.87(3) : 
3.56(2) : 1 for catenated rings. 

We investigate how the shape of rings changes when 
approaching the centers of mass, both for rings with- 
out topological interactions as well as for conformations, 
which obey the non-catenation constraint. Figure [9]A_ 
displays the results for a chain of length A^ = 1024. The 
polymers get strongly elongated for small ccnter-of-mass 
separations s — r/Rg. This elongation at full overlap 
s = remains in the limit of large chain length N ^ 00 
as is displayed in Fig. |9p. For non-catenated rings, the 
effect of elongation is slightly more pronounced than for 
rings, where topological effects are neglected. At full 
overlap we find for non-catenated rings that the ratios 
(A3) : (A2) : (Ai) are asymptotically 11.05(8) : 4.32(3) : 1, 
while for rings without mutual topological interactions 
we have 10.84(6) : 4.17(2) : 1. 

While the analysis of the radius of gyration indicated 
that the spatial extend of rings becomes larger, we have 
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0.02 



0.5 1 1.5 2 2.5 3 3.5 0.005 0.01 0.015 0.02 

CM separation r/Rgyr inverse chain length 1/N 

FIG. 9. A. Ratios between the gyration tensor's eigenvalues Ai < A2 < A3 in relation to the center of mass distance s = r/Rg. 
Data is shown for a system with chain length A'' = 1024 for two-ring conformations without mutual topological interactions 
(open circles) and non-catenated ring conformations (solid red circles). B. Finite-size behavior of the ratios (A3) / (Ai) and 
(A2) / (Ai) for a center-of-mass separation of s = 0. 



found here that this opening up does not lead to a com- 
plete mixing of the two rings, which would result in a 
more spherical structure; rather rings elongate strongly 
in one direction. 



F. Orientation and Alignment of rings 

To investigate the interactions between the two rings at 
full overlap in more detail, we look at the mutual align- 
ment and orientation. The alignment of the gyration 
ellipsoids of both rings with respect to each other can 
be investigated via two measures: Firstly, the average 
angle (cos 6*1) between the largest main axes of the el- 
lipsoids, secondly the average angle (cos6'2) between the 
vector connecting the centers of mass of both rings and 
the largest main axis of one ellipsoid. Its behavior has 
been observed for catenated rings^^^ \i -^^^g found that 
the main axes of the ellipsoids tend to align perpendicular 
for short CM separations, while for large CM separations, 
the alignment becomes more and more parallel. 

The results for two-ring conformations without mutual 
topological constraints and non-catenated ring conforma- 
tions as well as SAWs are shown in Fig.[TO] The black line 
at {cos9i) — 0.5 corresponds to the average angle for a 
random orientation of the vectors (Note that the angle 6i 
can adopt only values between zero and 90 degrees) . The 
average angle (cos^i) between the gyration tensors main 
axes (cf. Fig. 104.) displays a pronounced perpendicu- 
lar orientation for small CM separations r < Rg. This 



effect is markedly stronger for ring polymers than for 
linear ones. At intermediate separations of r « 1 — 2i?g, 
however, a regime is found where a slightly parallel align- 
ment of the rings is preferred compared to a random ori- 
entation. For large center of mass separations, a ran- 
dom orientation is approached; in this regime, the rings 
are nearly independent and therefore do not influence 



each other. An extrapolation to the infinite chain length 
regime at full overlap (s = 0) yields an average angle 
(cos^i) of 0.3543(7) for rings without topological inter- 
actions, 0.3492(6) for non-catenated rings and 0.3897(2) 
for self-avoiding walks. Details of the extrapolation are 
shown in the Supplementary Figure 3'*''. 



The angle (cos02) (see Fig. 10 3) shows a nearly ran- 
dom alignment for CM separations of r « 0. For interme- 
diate r « 1 — 2i?g, the alignment becomes pronouncedly 
more perpendicular compared to a random alignment. 
This effect is much stronger for non-catenated rings than 
for rings without mutual topological constraints. Clearly, 
ring polymers display a much stronger tendency to align 
non-randomly than self-avoiding walks. In the regime of 
perpendicular alignment at s — 1.2 the asymptotic val- 
ues of 0.3836(5) (rings without topological interactions), 
0.3661(4) (non-catenated rings) and 0.4327(3) (SAWs) 
are found for the average angle (cos 92^- 

A recurrent motif in the analysis of shape, dimensions 
and orientation is a change in structure from polymers 
at full overlap to polymers at intermediate separations 
(r ss 1 — 2i?g). While they are strongly elongated and 
aligned perpendicular at short separations, at intermedi- 
ate separations compaction and parallel alignment is ob- 
served. Clearly, this results from a tendency to minimize 
the overlap area between both rings. At full overlap, this 
is accomplished best by a strong elongation and perpen- 
dicular alignment (Fig. [Tl]^). When rings are separated 
further apart, the gyration ellipsoids can avoid intermin- 
gling by aligning in parallel (Fig. [TI)3). The restricted 
conformational space due to the presence of the other 
ring results in the observed compaction of the radius of 
gyration compared to isolated chains. 

Both the topological constraints and the ring connec- 
tivity of the chain have a strong influence on how the 
polymers are aligning when brought close together. In 
general, the view emerges that rings have a more aligned 
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0.5 1 1.5 2 2.5 3 
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FIG. 10. A. Average angle (cos Si) between the largest prin- 
cipal axes of the gyration ellipsoids of a two-chain conforma- 
tion in relation to its center of mass separations. The black 
line corresponds to a random orientation of the two axes. Re- 
sults are shown for self-avoiding walks (grey triangles), non- 
catenated rings (solid red circles) and rings without mutual 
interactions (open black circles). B. Average angle (cos 62) 
between the vector connecting the centers of mass and the 
largest principal axis of one polymer coil in relation of the 
CM separation. 



and ordered state than linear chains, an effect which is 
even amplified by the non-catenation constraint. 



IV. CONCLUSIONS 

In this study we investigated the effect of the non- 
catenation constraint on a system with two ring poly- 
mers. Our main focus was on the quantitative analysis 
of the strength of the interactions between ring polymers, 
including their topological interactions. For this pur- 
pose, we sampled conformations of isolated rings using 
a well-established Monte Carlo method. We evaluated 
the potential of mean force using the idea of effective 
interactions^^ following the method introduced by Daut- 
enhahn and HalP"^. The topological state of a two- ring 





center-of-mass 
distance s=1 .5 



FIG. 11. Sample conformations of non-catenated rings forced 
to a certain center-of-mass distance. The left-hand figure dis- 
plays a conformation of two non-catenated rings at nearly full 
overlap {s = r / Rg = 0.3). The two rings markedly align per- 
pendicularly. The right-hand figure shows a two-ring confor- 
mation of non-catenated rings at separation s = r/Rg = 1.5 
displaying a tendency for parallel alignment. 



conformation is analyzed by means of the Gauss linking 
number. 

We found that the effective potential at full overlap 
adopts a constant value both for linear chains and ring 
polymers (Figs.[2]and|3|. However, the repulsive interac- 
tion between the centers of mass is about 3 times larger 
for ring polymers than the effective potential for corre- 
sponding linear chains (Fig. [4| at small center-of-mass 
separations r. 

If the ring polymers have to stay in the fixed mutual 
topological state of non-catenation, the space of accessi- 
ble conformations is further reduced, thus the effective 
potential increases. We have evaluated the strength of 
these interactions - the topological potential Utop{r) - 
resulting purely from topological constraints. We find 
that the non-catenation constraint further increases the 
total repulsive interaction by about 10%. The strength 
of the potential is of the order of 1 /csT at small separa- 
tions r. The number of rejected conformations increases 
by over 50% compared to only taking into account ex- 
cluded volume interactions. 

The analysis of conformational properties of polymers 
brought close together reveals effects of the ring struc- 
ture both on size and shape. Ring polymers in proximity 
become swollen (Fig. |6| and strongly elongated (Fig. [9]) 
compared to the situation of isolated or far apart ones. 
While both effects are found for linear polymers, too, 
the effects are more pronounced for their ring counter- 
parts. Additional effects from the non-catenation con- 
straint concerning size and shape are visible, but small, 
indicating that the changes are mainly induced by the 
presence of additional material due to the rings being 
more compact than linear polymers. 

Remarkable effects of the ring structure are found con- 
cerning the alignment of ring polymers in proximity. 
There is a strong tendency to align perpendicular for 
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short center-of-mass separations, an effect which is much 
more pronounced than for hnear chains. At intermediate 
separations (r ^ 1 — IRg) a sHghtly parallel alignment is 
encountered. 

The transition from a linear to a ring polymer thus 
induces strong changes in the conformational and struc- 
tural properties of a system. In fact, ring polymers adopt 
a much more ordered and regular state, showing less in- 
termingling due to the increased repulsive interactions. 
The mutual alignment of ring polymers becomes much 
more different from a random state than found for linear 
polymers. 

Our findings suggest, that the ring topology of certain 
biopolymers like DNA and proteins, displays a benefit 
compared to the linear organization form. The forma- 
tion of chromatin loops in higher eukaryotes therefore 
might play a dominant role concerning overall nuclear or- 
ganizational principles. We have shown that loop struc- 
tures lead to a strong effective repulsion and for a first 
time quantitatively analyzed the strength of the result- 
ing interactions. Surely, the formation of multiple loops 
in the system of chrom atin which has been proposed in 
several modelJ^'^oIlIIMl jj^fj^ces even stronger effects and 
therefore might be responsible for maintaining the seg- 
regated state of chromosomes found in several experi- 
ment's]. Indeed, studies of looping polymers have re- 
vealed an eff ect on the abundance of inter-chromosomal 
contactJ^^HH. Whereas it has already been proposed 
that the ring structure prevents chromosomes from en- 
tangling^^ such a topology also induces a much more 
ordered and aligned state than linear polymers. It re- 
mains to be shown, to which extend similar results are 
found for chromatin models with loop^^flMI expect 
that both the local behaviour (e.g. inside of transcription 
factories or heterochromatic regions) and the global be- 
havior (segregation of chromosomes) might well be dom- 
inated by entropic and topological forces of chromatin 
loops. Contrary to a melt of ring polymers, which adopts 
a compact state, the dynamic formation of loop might 
provide a versatile mechanism for the cell to switch be- 
tween and maintain different local levels of compaction 
by facilitating or repressing loop formation. 
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